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SECTION  1 
INTRODUCTION 

SOFEPL  is  a  computer  program  for  drawing  line  plots  of  the 
output  from  SOFE,  a  Monte  Carlo  simulation  program.  SOFEPL 
(pronounced  so-fe-pe-el)  stands  for  Simulation  for  Optimal 
Filter  Evaluation  Plotter.  SOFEPL  is  flexible,  efficient,  and 
quickly  accessible.  This  user's  manual  covers  the  function, 
structure  and  operation  of  the  SOFEPL  program.  These  same  sub¬ 
jects  are  explained  for  SOFE  in  the  SOFE  user's  manual 
[Reference  1],  and  are  discussed  only  incidentally  here. 

Realizing  the  power  inherent  in  Monte  Carlo  simulation 
requires  an  extensive  statistical  computation  and  plotting  capa¬ 
bility.  Although  SOFE  offers  a  wide  range  of  options  for 
printing  results,  it  has  only  rudimentary  plotting  facilities  and 
no  statistical  computation  abilities.  These  plotting  and  sta¬ 
tistical  functions  are  accomplished  by  the  FORTRAN  program 
SOFEPL,  a  companion  postprocessor  to  SOFE.  The  interplay  of  SOFE 
and  SOFEPL  is  illustrated  in  Figure  1. 

This  user's  manual  describes  SOFEPL  functions,  documents 
the  program's  code  structure,  and  provides  instructions  for  and 
examples  of  SOFEPL  use.  Section  2  describes  the  problem  that 
SOFEPL  must  solve  and  gives  information  defining  its  abilities. 
Section  3  is  a  collection  of  program  notes  describing  such  topics 
as  program  organization,  data  structures  and  coding  conventions. 
Section  4  specifies  the  format  and  content  of  program  inputs  and 
outputs,  and  gives  instructions  for  program  use.  In  Section  5, 
representative  samples  of  SOFEPL ' s  capabilities  are  demonstrated 
for  the  orbit  problem  documented  in  [1].  Appendix  A  contains 
input  and  output  for  the  example  discussed  in  Section  5. 

Appendix  B  covers  DISSPLA  postprocessors.  (DISSPLA  is  the 


1 


Figure  1.  SOFE/SOFEPL  System  Diagram 


graphics  software  package  that  SOFEPL  uses.)  Appendix  C  descri¬ 
bes  two  merge  programs  that  combine  multiple  SOFE/SOFEPL  inter¬ 
face  files.  Readers  interested  in  exercising  the  program  in 
minimum  time  should  scan  Section  2,  skip  Section  3,  and  con¬ 
centrate  on  Sections  4  and  5. 


SECTION  2 

GENERAL  DESCRIPTION 


SOFE  is  a  general-purpose  Monte  Carlo  simulation  program 
that  facilitates  the  synthesis  of  a  Kalman  filter  for  a  multisen¬ 
sor  system,  and  aids  in  the  analysis  of  that  system's 
performance.  For  a  particular  system  configuration,  which  con¬ 
sists  of  both  truth  and  filter  models  and  measurement  updates 
from  various  sensors,  SOFE  executes  runs  over  a  specified  time 
interval  to  simulate  missions. 

As  each  simulation  run  develops,  SOFE  records  it  output  on 
an  unformatted  interface  file  designated  TAPE 4 .  This  file  accu¬ 
mulates  time-tagged  samples  of  the  five  vector  variables  named  in 
Table  1. 


TABLE  1 

VARIABLES  RECORDED  BY  SOFE  ON  TAPE4 


Description  of  Vector  Variable 

Symbol 

Dimension 

1 

Truth  state 

Xs 

NS 

2 

Filter  state  estimate 

Xf 

NF 

3 

Filter  variances 

Vf 

NF 

4 

Measurement  residuals 

ZRES 

M 

5 

Measurement  residual  variances 

ALPHA 

M 

Note  that 

Vf  is  the  diagnoal  of  the  filter 

covariance 

Pf.  Thus 

t  h 

each  element  Vf.  of  Vf  is  the  estimated  error  variance  of  the  i 

^  A 

filter  state  estimate  Xf^;  i.e.,  Vf^  =  Pf^. 


I 


The  data  of  Table  1  is  recorded  at  equally  spaced  intervals 
and  at  measurement  update  times.  This  sampling  scheme  permits 
the  user  to  capture  a  complete  history  of  each  Monte  Carlo  run. 


A  set  (ensemble)  of  runs,  taken  together,  comprises  a 
study.  The  statistical  computations  of  SOFEPL  are  accomplished 
at  each  sample  time  across  the  ensemble  of  runs  in  the  study. 
These  statistical  computations  consist  of  simple  averages  and 
standard  deviations  on  elements  (and  element  differences)  of  the 
five  variables  listed  in  Table  1.  SOFEPL  plots  these  averages 
and  standard  deviations  against  time  to  facilitate  analyzing  the 
performance  of  the  integrated  system. 


2.1  Average,  Standard  Deviation,  Estimation  Error 


Figure  2  illustrates  the  average  and  standard  deviation 
curves  that  SOFEPL  draws.  These  curves  derive  from  elementary 
statistical  computations.  To  avoid  misunderstandings  concerning 
these  computations,  they  are  developed  here. 


Consider  a  continuous  function  y(t)  which  can  have  a  dif¬ 
ferent  realization  on  each  (SOFE)  run.  Figure  2  illustrates  such 
a  function  for  three  runs,  y^t),  y2(t),  and  y3(t).  Of  course, 
y(t)  could  be  any  element  of  any  vector  in  Table  1.  Assume  that 
y(t)  is  sampled  on  each  run  at  time  t'.  If  there  are  a  total  of 
N  runs  in  the  study,  the  collection  Y(t')  of  N  samples  at  t'  is 


Y(t')  =  |yx(t ' ) ,y2(t' ) ,  ...  ,yN(t')| 


where  the  subscripts  n=l,2,...,N  represent  the  run  number.  In 
Figure  2,  N=3  runs.  The  arithmetic  mean  or  average,  y(t’),  of 
the  quantities  in  the  collection  Y(t')  is 


?(C>  -  i  l  yn(f) 
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Define  y(t^)  as  the  time-sequenced  collection  of  all  averages 
y(t')  for  all  sample  times  t'  in  the  closed  interval  [TO,TF] 
where  TO  and  TF  are  the  initial  and  final  time  of  each  run.  If 
*(t^)  is  shortened  to  y  for  notational  simplicity,  and  if  the 
sequence  of  sample  times  is  t^  ,  t2,  ...  ,  tR,  then  y  is 

y  =  |y(t1).  y(t2),  •••  'y^)}  (3) 

When  the  point  pairs  (t^,  y(t^))  of  (3)  are  computed,  plotted, 
and  connected  in  SOFEPL  by  straight  line  segments,  they  form  a 
curve  of  ensemble  averages  of  y(t)  samples.  The  dark  solid  line 
in  Figure  2,  connecting  K=4  sample  points,  is  such  an  average. 

In  all  that  follows,  the  overbar  notation  will  represent  a  time 
history  as  in  (3)  of  averages  formed  over  the  runs  in  the  study, 
each  run  weighted  equally  as  in  (2). 

SOFEPL  also  computes  a  standard  deviation,  denoted  here  by 
a  capital  S.  For  the  set  of  N  samples  in  (1),  the  standard 
deviation  is 


where 


S(f) 


4T 


&t]V 


V(t’  ) 


V(t’ ) 


N 

N-l 


y2 ( t  * ) 


(4) 


(5) 


Here  V(t')  is  the  familiar  variance  computation.  Equation  (4)  is 
an  approximate  result  that  appears  in  Reference  2,  page  76.  A 
full  derivation  of  the  exact  result  appears  in  Reference  3, 
Chapter  15.  The  advantage  of  (4)  is  that  it  produces  an  unbiased 
estimate  of  S(t’),  especially  when  N^9  and/or  when  y(t')  is  a 
sample  from  a  normal  or  nearly  normal  distribution.  These  latter 
two  conditions,  however,  need  not  be  rigorously  satisfied. 

The  discrete-time  function  S(t^),  or  S  for  short,  is  the 
time-sequenced  collection  of  all  standard  deviations  S(t')  for  t' 
in  [TO, TF] .  Thus,  in  analogy  with  (3), 
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The  dash-dot  line  in  Figure  2  is  the  S  curve  of  y(t)  for  the 
three  runs  shown.  Note  that  S  is  not/ defined  for  N=1  run. 

Define  E.  .  as  the  estimation  error  for  the  jth  state  in  the 
*  J  A 

estimated  state  vector  Xf.  Then  E^  is 


(7) 


where  Xs^  and  Xf  ^  are  the  1  and  j  elements  of  Xs  and  Xf, 
respectively.  In  (7),  Xs.  is  chosen  to  force  the  comparison  with 

A  * 

Xfj  to  be  between  the  same  physical  quantities.  SOFEPL  provides 
a  means  for  forming  any  difference  j ,  l<i<NS,  l<j<NF. 


2 . 2  Available  Plots 

Table  2  lists  the  18  plots  that  SOFEPL  can  draw.  All  plots 
are  time  histories  of  the  variables  shown  in  the  table,  all 
averages  are  ensemble  averages  per  (3),  and  all  standard 
deviations  are  formed  per  (6).  A  representative  set  of  plots 
from  Table  2  is  described  in  the  following  paragraphs. 

j.  u 

In  plot  type  1,  the  i  state  of  the  truth  state  vector  Xs 
is  averaged  over  the  ensemble  of  runs  in  the  study  and  plotted 
through  time  to  produce  a  single  curve.  In  type  3,  the  average 
of  the  ith  state  of  Xs  forms  one  curve  while  the  average  of  the 
jth  state  of  Xf  forms  the  other. 

In  type  6,  the  average  of  the  difference  of  truth  state 
i  and  filter  state  j  (see  (7))  forms  one  curve;  the  standard 
deviation  SE^  of  this  difference,  as  prescribed  by  (6),  is  the 
second  curve,  and  -SE,  .  is  the  third  curve.  A  type  9  plot  is  the 

^  f  h 

average  of  the  square  root  of  the  i  diagonal  element  Vf^  of  the 
filter  covariance  Pf;  this  is  the  filter's  uncertainty  about  its 
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TABLE  2 

AVAILABLE  PLOT  TYPES 

Plot 

Type 

No.  of 
Curves 

Description  of 

Curves 

1 

1 

Xs. 

2 

1 

A 

Xf. 

3 

2 

*si 

If. 

3 

4 

3 

A  _ 

Xf  . 

1 

*fi+/v£7 

5 

1 

E.  . 

13 

6 

3 

E.  . 

13 

+SE .  . 

13 

-SE.  . 

13 

7 

3 

'Eij 

E.  ,+SE.  . 

13  13 

E. .-SE. . 

13  13 

8 

3 

~Eij 

+/VF7 

3 

-/ vf . 

3 

9 

1 

/VF7 

l 

10 

1 

SE.  . 

11 

1 

12 

2 

SE.  . 

13 

/w- 

13 

2 

/sEij+Eij2 

✓  V£  . 

3 

14 

5 

Lj 

E.  .  +  SE.  . 
13  -  13 

+/VFT 

15 

1 

ZRES. 

1 

16 

3 

ZRES . 

i 

+✓ ALPHA. 
i 

-/ALPHA'" 

i 

17 

2 

sz . 

1 

/alpha: 

i 

18 

NCRV 

User  defined 

using  arrays 

P1-P10 

estimate  Xfi.  Note  that  standard  deviations  Vvf . ,  not  variances 
Vf .  ,  are  averaged.  Plot  type  11  is  the  root  mean  square  (RMS) 

t  A 

value  of  the  estimation  error  associated  with  filter  state  Xf  . . 

3 

In  type  16,  the  average  of  the  ifc^  of  M  residuals  forms  one 
curve.  The  remaining  two  curves  are  plus  and  minus  the  average 
of  the  square  root  of  the  ifc^  element  in  the  vector  ALPHA. 

ALPHA ^  is  the  filter's  predicted  uncertainty  for  measurement 
residual  ZRES^.  Plot  type  17  compares  this  predicted  value  to 
SZ^ ,  the  actual  standard  deviation  of  the  ZRES^  data  over  the  N 
runs  in  the  study.  SZ^  is  computed  using  (6),  as  was  SE^j. 

Plot  type  18  is  a  user-defined  plot.  Whereas  types  1 
through  17  are  entirely  determined  by  the  parameters  in  the 
user's  input  data,  type  18  is  formed  in  large  part  by  FORTRAN 
code  constructed  by  the  user.  Type  18  is  discussed  in  detail  in 
Subsection  4.2. 

2 . 3  Other  Attributes 

SOFEPL  is  written  in  FORTRAN  and  was  developed  on  Control 
Data  Corporation  (CDC)  equipment  where  it  loads  in  approximately 
115,000  octal  words  of  memory.  A  concerted  effort  has  been  made 
to  use  programming  conventions,  1966  ANSI  constructs,  modular 
concepts,  and  informative  comments  that  will  expedite  the 
transfer  of  the  program  to  other  machines.  In  addition,  many 
measures  have  been  taken  to  make  the  program  efficient  in  its  use 
of  central  memory  and  execution  time. 

SOFEPL  calls  subroutines  in  the  graphics  package  DISSPLA 
for  the  actual  plot  generation.  The  output  of  this  package  is  a 
compact  file  of  plotting  commands  in  a  special  format.  This 
file,  called  PLFILE,  is  then  submitted  to  a  DISSPLA  postprocessor 
program  to  make  the  actual  plots  on  the  desired  plotting  device. 
While  it  may  seem  redundant  for  a  postprocessor  (SOFEPL)  to 
require  a  second  postprocessor  to  complete  the  final  product, 
this  approach  separates  plot  file  generation  from  plotting,  and 
thereby  provides  flexibility. 
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The  following  peripherals  are  required  to  run  SOFEPL: 

•  Card  reader  or  data  input  terminal  to  enter  input  para¬ 
meters. 

e  A  tape  or  disk  file  to  supply  TAPE4  data  from  SOFE. 

•  A  line  printer  to  print  listable  output. 

•  A  tape  or  disk  file  to  store  plot  commands  (PLFILE)  for 
subsequent  plotting  by  the  DISSPLA  postprocessor. 

DISSPLA  postprocessors  in  turn  usually  require  a  similar  set  of 
peripherals  for  reading  the  plot  commands  and  writing  to  the 
plotting  device. 


SECTION  3 
PROGRAM  NOTES 

SOFEPL  is  a  tool  for  plotting  results  from  a  series  of  SOFE 
runs.  From  a  user's  point  of  view,  the  program  is  flexible  in 
the  kind  of  plots  it  can  make,  in  the  choice  of  runs  it  uses,  and 
in  its  presentation  format.  In  addition,  SOFEPL  executes  quickly 
and  uses  core  efficiently.  This  section  describes  the  code 
structure  and  some  key  features  of  the  program  design  that 
contribute  to  accomplishing  the  goals  listed  above. 

3 . 1  Organization 

SOFEPL  is  organized  into  five  major  sections,  each  section 
representing  a  phase  in  the  process  of  creating  a  plot(s)  from 
SOFE  output.  These  five  sections  are  designated  as  follows: 

(1)  Program  initialization, 

(2)  Plot  selection, 

(3)  Data  accumulation, 

(4)  Plot  execution, 

(5)  Program  completion. 

The  organization  of  these  five  tasks  is  shown  in  the  macro-level 
flowchart,  Figure  3. 

Program  initialization  creates  data  and  storage  areas 
inside  SOFEPL  based  on  information  from  the  header  of  the  inter¬ 
face  file,  TAPE4.  Plot  selection  reads  the  user's  input  cards 
which  define  each  plot  desired  and  sets  aside  the  necessary 
storage  areas  in  memory.  The  data  accumulation  section  reads 
data  from  TAPE4  and  stores  that  data  in  the  appro*  riate  storage 
areas.  Plot  execution  calculates  the  desired  curves  from  the 
stored  data  and  constructs  a  file  containing  the  requested  plots. 
Program  completion  closes  the  plot  file  and  ends  the  run. 

SOFEPL  accomodates  plots  with  varying  numbers  of  curves  and 
points  by  employing  a  'spooling'  technique  for  data  storage. 
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With  this  technique,  the  amount  of  storage  allocated  to  each  plot 
is  exactly  what  is  needed  for  that  plot  and  waste  of  core  storage 
is  held  to  a  minimum. 

SOFEPL's  spooling  technique  works  as  follows.  The  spool  is 
a  single  linear  array  named  ‘A*  of  dimension  10000.  As  each  plot 
request  is  read,  it  is  interpreted  (in  light  of  information  from 
TAPE4 1 s  header)  to  determine  this  plot's  storage  requirements. 
Blocks  of  storage  are  then  set  aside  in  A  to  hold  the  points  for 
each  curve  as  well  as  other  necessary  data.  These  blocks  are 
accessed  using  an  index  of  pointers.  Furthermore,  the  block  for 
plot  n  is  adjacent  to  that  for  n-1  and  n+1  with  no  unused  words 
between  them;  thus  loss  of  core  storage  is  minimized. 

When  more  than  one  plot  is  requested,  execution  time  is 
also  reduced  by  use  of  the  spooling  technique.  This  occurs 
because  several  plots  are  being  made  on  each  pass  through  TAPE4, 
thereby  reducing  the  number  of  times  TAPE4  is  read  and  saving  I/O 
time.  Clearly,  if  the  spool  was  long  enough,  all  plots  could  be 
made  in  one  pass  through  TAPE4 .  The  tradeoff  is  between  the  cost 
of  core  storage  and  the  cost  of  I/O  time.  This  tradeoff  is 
approximately  correct  for  the  CDC  system  with  A  dimensioned 
10000.  This  dimension  is  easily  changed  for  another  system,  or 
for  plots  with  huge  numbers  of  points,  by  altering  a  COMMON  card 
and  a  DATA  card  in  SOFEPL's  main  routine. 

3 . 2  Code  Structure 

SOFEPL  consists  of  a  main  executive,  18  subroutine 
subprograms,  and  two  function  subprograms.  Additionally,  SOFEPL 
references  nine  DISSPLA  subroutines,  seven  FORTRAN  library  func¬ 
tions,  and  three  system  utility  functions.  The  main  executive 
sequences  the  program  through  its  five  sections  by  calling  on  the 
appropriate  routines  for  each  plot-creation  phase. 

Figure-  4  indicates  the  flow  of  SOFEPL  through 
its  routines.  This  figure  illustrates  which  routines  accomplish 
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Figure  4.  Subroutine  Calling  Sequence. 
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which  solution  phase  and  defines  dependencies  in  the  collection 
of  routines.  None  of  the  FORTRAN  library  functions  that  SOFEPL 
uses  are  shown  in  Figure  4.  A  subprogram  name  prefixed  by  an 
asterisk  indicates  a  DISSPLA  subroutine  or  system  utility 
function,  while  the  remaining  names  are  SOFEPL  routines. 

The  first  section,  program  initialization,  begins  with  a 
call  to  INITAL  where  the  TAPE4  header  is  read  to  acquire  data 
that  will  allow  arrays  to  be  sized  properly.  INITAL  then  writes 
the  SOFEPL  header  page  and  calls  SPLITA  where  ‘offset  pointers' 
are  assigned  in  A  to  locate  quantities  and  data  blocks.  INITAL 
checks  to  be  sure  that  the  spool  size  of  10000  words  is  large 
enough  to  handle  the  user's  problem;  if  not,  a  warning  message  is 
written  and  the  program  is  STOPed.  Finally,  INITAL  initializes 
the  DISSPLA  software  by  calling  COMPRS. 

During  the  second  phase,  plot  selection,  subroutine  CARDIN 
takes  over  and  calls  on  CARD12  to  read  cards  1  and  2,  and  on 
CARD35  to  read  cards  3,  4,  and  5.  Subroutine  CHKDTA  is  called  by 
CARD12  to  range-check  the  parameter  values  read  from  card  1,  the 
plot  selection  card.  Subroutine  STRING  is  called  by  CARD12  to 
acquire  the  desired  runs  from  card  2,  the  run  selection  v_  rd. 
CARD35  reads  the  label  data  for  the  plot  and  [for  the  y-  and 
t-axes.  The  routines  in  this  section  are  '  v  -atec  -.ib'hil  the 
spool  is  filled  with  plots  or  until  all  p^.t  sets  are  read, 
whichever  comes  first. 

When,  in  the  data  accumulation  phase,  a  new  run  is  ready 
for  processing,  subroutines  PLTNED  and  SCAN  are  called  to  deter¬ 
mine  which  plots  require  the  data  from  this  run.  If  no  plots 
require  this  run's  data,  SKPRUN  is  called  to  skip  over  the  run 
and  by-pass  its  data  entirely.  Otherwise  STORE  is  used  in  con¬ 
junction  with  USTORE  and  P15T17  to  read  and  accumulate  the  data 
in  the  correct  array  locations. 

The  fourth  section,  plot  execution,  uses  SETPLT  in  conjunc¬ 
tion  with  USETPL  to  form  averages  and  standard  deviations  from 
the  accumulated  data  arrays.  PLOT  is  then  called  to  draw  the 
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requested  plots.  Seven  of  the  nine  calls  to  DISSPLA  routines 
occur  in  PLOT. 

DONEPL  is  called  as  the  final  act  of  SOFEPL  execution  to 
close  PLFILE .  This  last  phase  and  the  first  phase  occur  only 
once  for  each  SOFEPL  run,  while  the  middle  three  phases  iterate 
until  the  entire  plot  specification  file  is  processed. 

3 . 3  Coding  Conventions 

SOFEPL  is  written  in  FORTRAN  using  primarily  1966  ANSI 
standard  constructs.  Exceptions  include  use  of  quoted  comments 
after  STOPs,  arithmetic  expressions  in  array  arguments,  list 
directed  input,  IMPLICIT  statements,  quoted  hollerith  data  in 
FORMATS,  and  a  few  others  common  to  most  FORTRAN  compilers. 

In  writing  SOFEPL  we  attempted  to  use  clear,  easily- 
understood  code  that  is  organized  in  top-down  fashion.  Nearly 
40  percent  of  the  statements  are  comments  that  are  designed  to  be 
concise,  informative,  and  accurate.  SOFEPL  uses  only  single  pre¬ 
cision  variables  of  the  REAL,  INTEGER,  and  LOGICAL  type. 

Variables  are  given  meaningful  names  from  the  36  alphanumeric 
characters .  The  first  letter  of  a  variable ’ s  name  conforms  to 
the  FORTRAN  default  rule  for  integer  and  real  names.  (STRING, 
STRGET,  and  STRNTP  are  exceptions  where  all  variables  are  forced 
to  be  integers  through  IMPLICIT  statements.)  Counters  generally 
start  with  I  or  J,  pointers  with  K,  logical  variables  and  labels 
with  L,  and  sizes  and  dimensions  with  M  or  N.  All  format  state¬ 
ments  are  located  at  the  end  of  routines. 

3 . 4  SOFEPL  Data  Structures 

SOFEPL  uses  two  major  data  structures.  The  first  is  the 
SOFE/SOFEPL  interface  file,  TAPE4.  The  other  is  the  A/lA  array 
used  to  store  the  required  data  for  the  selected  plots. 


18 


3.4.1  SOFE/SOFEPL  Interface  File 

The  SOFE/SOFEPL  interface  file  is  an  unformatted 
sequential  file  written  in  internal  binary  representation  by 
SOFE.  This  file,  called  TAPE4,  contains  a  header  record  followed 
by  IPASS  subfiles,  where  IPASS  is  the  number  of  SOFE  simulation 
runs.  Each  subfile  (referred  to  in  subsequent  discussion  as  just 
a  "file")  is  separated  from  its  neighbor  by  an  end-of-file  mark. 
Figure  5  shows  the  overall  structure  of  TAPE4. 

The  header  record  of  TAPE4  contains  data  pertaining 
to  the  originating  SOFE  run,  namely  the  values  of  NS,  NF,  M,  TO, 
TF,  TMEASO,  DTMEAS,  DTCCPL,  IPASS,  run  date,  run  time,  and  run 
title.  This  record  is  read  and  echoed  in  INITAL.  Each  file  on 
the  interface  tape  contains  a  series  of  variable  length  records. 
In  turn  each  record  is  comprised  of  vector  variables  sampled  at 
an  event  time  as  specified  in  Table  3.  The  records  for  file  one 
are  output  during  SOFE  run  one,  those  from  file  two  during  run 
two,  and  so  on.  The  event  times,  which  are  determined  by  SOFE 
parameter  values,  are  identical  in  every  run.  Records  are  read 
and  processed  in  subroutine  STORE  where  the  read  process  is  aided 
by  the  fact  that  the  record's  length  is  its  first  word.  Note  in 
Table  3  that  t^  may  equal  t^.  When  this  occurs,  a  double  sample 
of  Xs,  Xf  and  Vf  occurs  at  t^  thereby  causing  these  points  to  be 
double  plotted. 

3.4.2  The  A/lA  Array 

Because  the  A  array  is  equivalenced  to  the  IA  array 
so  that  both  integers  and  real  values  are  readily  accessible, 
the  name  is  A/IA.  A/lA  was  referred  to  earlier  as  the  spool  or 
spooling  area.  Most  of  the  variables  and  arrays  necessary  for 
generation  of  plots  are  stored  in  this  array.  These  variables 
and  arrays  are  referenced  by  means  of  an  index  of  pointers. 
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Figure  5.  TAPE 4  Structure. 
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TABLE  3 

TIMING  OF  RECORDING  OF  VECTOR  VARIABLES  ON  TAPE4 


Vectors 

Time  of  Sample 

Number  of 

Recorded 

Symbol 

Description 

Samples 

Per  Run 

Xs  Xf  Vf 

- 

TO 

Initial  time 

1 

Xs  Xf  Vf 

Fixed  sampling  interval 
( =DTCCPL ) 

(TF-TO) 

"  DTCCPL 

Xs  Xf  Vf 

X 

Before  measurement 
incorporation 

(TF-TMEASO) 

DTMEAS 

A 

Xs  Xf  Vf 

+ 

t. 

1 

After  measurement 
incorporation 

II 

Xs  Xf  Vf 
ZRES  ALPHA 

t.+c 

1 

After  measurement 
incorporation  and 
feedback  correction 

II 

A 

Xs  Xf  Vf 

TF 

Final  time 

1 

Figure  6  shows  the  general  structure  of  the  A/IA  array.  The 
number  of  plots  for  which  data  can  be  stored  is  determined  by  the 
dimension  of  A/IA,  by  the  curves  per  plot,  and  by  the  number  of 
points  per  curve.  If  there  is  more  storage  available  in  A/lA 
than  necessary,  then  only  one  pass  through  the  interface  file  is 
made.  If  more  plots  are  requested  than  can  be  accomodated  in  the 
A/lA  array  on  one  pass,  then  plot  data  for  as  many  plots  as  will 
fit  is  stored.  When  plot  set  reading  is  complete,  the  program 
makes  a  pass  through  the  interface  file  to  generate  the  actual 
plots.  When  the  first  set  of  plots  is  finished,  the  data  in  A/lA 
for  these  plots  is  'dumped'  so  SOFEPL  can  reuse  this  storage 
space  for  the  next  set  of  plots,  and  so  on  until  all  plot 
requests  are  processed. 
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The  first  section  of  the  A/IA  array  is  the  general 
work  area.  Figure  6  shows  the  internal  structure  of  the  general 
work  area  and  its  relationship  to  the  rest  of  the  array. 

Figure  7  shows  the  structure  within  a  plot  data  area.  It  should 
be  noted  in  Figure  7  that  the  arrays  Y1-Y5  and  P1-P10  are  allo¬ 
cated  storage  only  for  the  plot  types  for  which  they  are  needed. 
Storage  is  allocated  by  adjusting  the  array  pointers  in 
subroutine  CARDIN.  For  example,  suppose  a  plot  type  4  is 
requested,  requiring  three  Y  arrays,  Yl,  Y2,  Y3 .  Suppose  that 
after  allocating  storage  for  the  preceding  variables,  the  storage 
for  the  Y  arrays  starts  at  A/lA(1000).  This  means  that  the  last 
value  in  ISETS1/ISETS2  is  at  location  A/lA(999).  Also  suppose 
that  a  run  on  the  interface  tape  contains  150  data  points.  Since 
the  storage  for  the  Y  arrays  begins  at  A/lA(1000),  SOFEPL  would 
place  1000  in  the  variable  KYI,  1150  in  KY2  and  1300  in  KY3. 

KYI,  KY2,  and  KY3  now  point  to  the  locations  of  Yl ,  Y2,  and  Y3  in 
the  A/lA  array.  Values  of  1300  would  also  be  placed  in  KY4,  KY5 , 
KPl ,  KP2 , . . . , KP10 ,  thus  implying  that  these  arrays  contain  zero 
storage.  If,  for  example,  SOFEPL  needed  to  access  the  tenth 
point  in  Y2,  the  following  code  would  be  used: 

(value  of  point  #10  in  Y2 )  =  A(lA(KLOCAT  +  KY2 )  +10-1) 

where 

KLOCAT  =  KPLPTR  +  IPLOT  -  1. 

The  complexity  of  the  pointers  is  the  price  paid  for  program  ver¬ 
satility  and  efficiency. 
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SECTION  4 

MAKING  PLOTS  WITH  SOFEPL 


The  procedure  for  obtaining  SOFEPL  plots  is  illustrated  in 
Figure  1  and  summarized  here: 

(1)  Run  SOFE  with  control  parameter  LCC  made  true  and  save 
TAPE4 ,  the  SOFE/ SOFEPL  interface  file. 

(2)  Construct  a  SOFEPL  'plot  specification  file*. 

(3)  Write  user  routines  USTORE  and  USETPL  if  a  type  18 
plot  is  desired. 

(4)  Run  SOFEPL  using  TAPE4  and  the  plot  specification  file 
as  inputs.  Save  the  output  file  PLFILE. 

(5)  Make  plots  on  a  plotting  device  using  PLFILE  as  input 
to  the  DISSPLA  postprocessor  for  that  device. 

Step  1,  the  process  of  running  SOFE  and  saving  TAPE4,  is 
discussed  thoroughly  in  [1].  Since  TAPE4  can  be  a  very  large 
file,  the  user  must  exercise  good  judgement  in  his/her  choice  of 
media  for  saving  it.  Magnetic  tape  is  usually  least  expensive 
while  disk  is  most  expensive.  Steps  2,  3,  and  4  of  the  above 
procedure  are  discussed  next  in  Subsections  4.1,  4.2,  and  4.3. 
Step  5  instructions  are  given  in  Appendix  B. 

4 . 1  The  SOFEPL  Plot  Specification  File 

The  SOFEPL  plot  specification  file  determines  which  plots 
are  produced  by  SOFEPL.  This  file  takes  the  form  of  user- 
supplied  data  cards,  or  card  images.  The  plot  specification  file 
consists  of  a  run  title  card  followed  by  any  number  of  five-card 
groups  (called  plot  sets)  which  specify  each  particular  plot. 

That  is,  SOFEPL  produces  one  plot  for  each  plot  set  submitted  by 
the  user  through  the  plot  specification  file.  Figure  8  shows  the 
overall  structure  of  the  plot  specification  file.  As  this  figure 


shows,  the  plot  specification  file  contains  these  six  types  of 
data  cards: 


T 


•  Run  Title  Card 

•  Plot  Specification  Card 

•  Run  Selection  Card 

©  Plot  Title  Card 

•  Time  Axis  Label  Card 

•  Y-Axis  Label  Card. 

The  next  four  subsections  describe  these  cards  in  detail. 

4.1.1  SOFEPL  Run  Title  Card 

The  first  card  of  the  plot  specification  file  is 
the  run  title  card.  This  card  is  used  to  input  a  title  that 
characterizes  the  upcoming  SOFEPL  run.  This  title  appears  only 
on  the  printed  output  provided  by  SOFEPL.  Thus  this  title  can  be 
used  to  match  the  printed  output  with  a  plot  specification  file. 
The  run  title  card  can  contain  up  to  80  characters.  This  card 
has  no  format  or  character  type  restrictions. 

4.1.2  Plot  Specification  Card 

The  first  card  of  a  plot  set  is  the  plot  specifica¬ 
tion  card.  This  card  is  used  to  input  the  15  parameters  that 
define  the  particular  plot  to  be  drawn.  These  parameters  are 
listed  in  Table  4;  their  functions  will  be  illustrated  in 
Section  5. 

The  parameter  values  on  this  card  are  read  by 
SOFEPL  using  CDC  list  directed  format.  Values  input  under  this 
CDC  free-form  convention  must  be  in  proper  order  (that  shown  in 
Table  4)  but  need  not  reside  in  preassigned  columns  on  the  card. 
In  addition,  multiple  quantities  may  be  entered  on  a  single  card 
separated  only  by  commas  or  blanks.  To  illustrate,  one  could 
enter  values  of  14,1,3,3  for  IPLTYP,  NXS,  NXF,  NVF  in  any  of 
these  ways: 

14,1,3,3/  or  14  1  3  3/  or 

14. . 2*3/  or  14,  ,3,3/  or 

14. .  3  3/  or  ... 
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TABLE  4 


PLOT  SPECIFICATION  CARD  PARAMETERS 


Position  Variable  Variable  Default  Data 


on  Card 

Name 

Def inition 

Limits 

Value 

Type 

1 

IPLTYP 

Plot  Type  (Table  2) 

1-18 

none 

Integer 

2 

NXS 

Truth  State  Element 
Number 

1-NS 

1 

Integer 

3 

NXF 

Filter  State  Element 
Number 

1-NF 

1 

Integer 

4 

NVF 

Covariance  Matrix 
Diagonal  Element 
Number 

1-NF 

1 

Integer 

5 

NZRES 

Measurement  Residual 
Number 

1-M 

max (1, NALPHA) 

Integer 

6 

NALPHA 

Measurement  Residual 
Variance  Number 

1-M 

NZRES 

Integer 

7 

IPNTSW 

Print  Switch 
(1  *  print  on) 

0  or  1 

0 

Integer 

8 

SFT 

T-Axis  Scale  Factor 

none 

1.0 

Real 

9 

SFY 

Y-Axis  Scale  Factor 

none 

1.0 

Real 

10 

TB 

T-Axis  Min.  Value 

none 

min  t 

Real 

11 

TE 

T-Axis  Max.  Value 

none 

max  t 

Real 

12 

YLO 

Y-Axis  Min.  Value 

none 

min  y 

Real 

13 

YHI 

Y-Axis  Max.  Value 

none 

max  y 

Real 

14 

NCURV 

Number  of  Curves  to 
Be  Drawn  (Plot  Type 
18  Only) 

1-5 

0 

Integer 

15 

NPS 

Number  of  P  Arrays 
Used  (Plot  Type  18 

NCRV-10 

0 

Integer 

Only) 
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Note  that  a  pair  of  commas  separated  by  an  optional  blank (s)  has 
the  effect  of  skipping  the  reading  of  the  parameter  at  that 
location;  thus  this  parameter  is  assigned  its  default  value  which 
is  1  for  NXS  (see  Table  4).  Also  note  the  repetition  factor  in 
the  third  version.  The  slash  at  the  end  of  each  version 
instructs  the  program  to  discontinue  reading  the  plot  specifica¬ 
tion  card;  the  parameters  remaining  in  the  input  list  then 
acquire  their  default  values.  Real  numbers  are  entered  using  any 
valid  FORTRAN  format.  The  user  can  consult  his/her  FORTRAN 
manual  for  additional  list  directed  information  but  the  examples 
above,  coupled  with  those  for  the  sample  in  Section  5,  should  be 
adequate  for  most  purposes. 

A  few  things  should  be  mentioned  about  the  scale  factor  and 
axis  limit  parameters  (SFT,  SFY,  TB,  TE,  YLO,  YHI ) .  When  time 
limits  (TB  and/or  TE)  are  specified,  the  points  lying  outside 
the  specified  range  do  not  appear  on  the  plot.  These  points  are, 
however,  transfered  to  user-defined  subroutine  USTORE  if  a 
type  18  plot  is  in  process.  The  same  is  true  for  plots  with  y- 
axis  limits.  All  points  are  used  in  the  various  SOFEPL 
calculations.  The  out-of-bounds  points  are  deleted  by  the 
DISSPLA  plotting  routines  which  are  the  last  called.  A  list  of 
the  out-of-bounds  points  is  generated  in  SOFEPL  printed  output. 
When  entering  both  a  scale  factor  and  axis,  limits,  the  axis 
limits  should  be  in  the  units  of  the  rescaled  data. 

4.1.3  The  Run  Selection  Card 

The  second  card  of  a  plot  set  is  the  run  selection 
card.  This  card  supplies  a  list  of  the  SOFE  simulation  runs  that 
are  to  be  used  in  forming  averages  and  standard  deviations  for 
this  plot.  This  card  contains  a  list  of  data  items  separated  by 
commas.  These  items  can  either  be  integers  reperesenting  par¬ 
ticular  runs,  or  integer  pairs,  separated  by  a  dash,  repre¬ 
senting  an  inclusive  set  of  SOFE  runs.  The  letter  'L'  can  also 
be  used  in  place  of  the  integer  value  of  the  last  SOFE  run.  For 
example,  if  there  were  17  runs  on  TAPE4,  a  card  reading  1,  5-9, 
11,  12,  15-L  would  instruct  SOFEPL  to  use  SOFE  runs  1,  5,  6,  7, 
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8,  9,  11,  12,  15,  16,  and  17.  In  this  case  N,  the  number  of  runs 

in  the  study,  is  11.  Note  that  the  list  of  data  items  must  be 
contained  on  one  card  and  must  increase  positively  across  the 
card.  Blanks  on  this  card  are  ignored. 

4.1.4  Plot,  Time-Axis,  and  Y-Axis  Title  Cards 

The  final  three  cards  of  the  five-card  plot  set  are 
the  plot  title,  the  time-axis  label,  and  the  y-axis  label.  Like 
the  run  title  card,  these  cards  are  read  with  a  standard  FORTRAN 
formatted  read.  The  titles  on  these  cards  may  be  up  to  60 
characters  long  and  should  be  composed  from  the  standard  alpha¬ 
numeric  characters,  A-Z  and  1-9,  plus  the  following  nine  special 
characters:  ',.+-*/().  A  $  placed  after  any  title  will 

cause  the  title  to  be  centered  in  the  appropriate  space. 

4 . 2  Using  Plot  Type  18 

Plot  type  18  enables  the  user  to  construct  new  variables 
from  the  elements  of  Xs,  Xf  and  Vf,  and  then  to  draw  time 
histories  of  up  to  five  such  variables  on  a  single  plot.  In  the 
example  to  be  presented  in  Section  5,  polar  coordinates  are  con¬ 
verted  to  rectangular  coordinates  through  the  sine  and  cosine 
relationships.  The  range  of  possibilities  for  plot  type  18  is 
obviously  very  broad.  It  is  hoped  that  this  capability  will 
satisfy  most  of  the  unusual  needs  not  covered  by  types  1-17. 

The  user  defines  the  variables  that  SOFEPL  will  plot  by 
writing  subroutines  USTORE  and  USETPL.  These  routines  are  only 
stubs  in  SOFEPL' s  code  until  the  user  inserts  executable  code  in 
them.  USTORE  operates  as  the  variable  computation  routine.  Each 

A 

time  new  samples  of  Xs ,  Xf,  and  Vf  are  read  from  TAPE4,  USTORE  is 
called  with  those  three  vector  samples  as  formal  parameters.  K 
such  calls  to  USTORE  occur  for  each  selected  run  and  N*K  calls 
occur  for  a  study.  (Note  that  K  and  N  are  denoted  KREC  and 
RUNCNT  in  the  stubs  for  USTORE  and  USETPL. )  In  USTORE,  the  user 
computes  the  variables  he  needs  (up  to  ten)  and  saves  them  in 
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appropriate  arrays.  Often  he  will  be  saving  sums  and  sums- 
squared  of  these  variables  in  order  to  form  averages  and  standard 
deviations  in  the  final  plots. 

After  all  N  runs  have  been  processed  through  USTORE,  SOFEPL 
calls  USETPL ,  thereby  affording  the  user  an  opportunity  for  some 
final  computations  and/or  rearrangements  before  plotting  begins. 
As  stated  above,  this  is  where  averages  and  standard  deviations 
should  be  computed  by  applying  (2)  and  (4)  at  each  sample  time. 
Since  all  of  SOFEPL' s  abilities  for  axis  scaling  and  windowing 
remain  in  effect  for  plot  type  18,  USETPL  is  unburdened  from 
these  duties. 

SOFEPL  provides  the  type  18  user  with  up  to  ten  scratch 
arrays  (PI,  P2,  ...  ,  P10).  The  actual  number  of  arrays  provided 

is  NPS,  an  integer  entered  on  the  plot  specification  card  espe¬ 
cially  for  plot  type  18.  These  arrays  are  dimensioned  K  long, 
large  enough  to  hold  one  value  for  every  sample  time  in  a  SOFE 
run.  Curves  can  be  drawn  from  the  first  five  of  these  P  arrays. 
The  actual  number  of  curves  drawn  is  NCRV,  an  integer  also 
entered  on  the  plot  specification  card,  especially  for  this  plot 
type.  Additional  storage  arrays  can  be  obtained  via  additional 
DIMENSION  statements  within  the  USTORE  and  USETPL  routines. 

The  definitions  for  the  formal  parameters  in  the  argument 
lists  of  USTORE  and  USETPL  are  given  in  Tables  5  and  6.  It 
should  be  noted  that  all  ten  P  arrays  are  named  in  the  argument 
lists  of  both  USTORE  and  USETPL  even  though  the  user  has 
requested  only  NPS  of  them  (NPS  j<  10)  .  NPS  is  used  to  determine 
which  of  the  P  arrays  are  allocated  storage.  Arrays  PI,  P2, 

...  ,  PNPS  are  allocated  storage,  whereas  the  arrays  P(NPS+1), 
P(NPS+2),  ...  ,  P10  are  not.  Thus  this  later  set  of  P  arrays 

must  not  be  used. 

4 . 3  JCL  for  Running  SOFEPL 

The  job  control  language  (JCL)  constructs  required  to  run 
SOFEPL  on  the  CDC  system  at  Wr ight-Patterson  Air  Force  Base  will 
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TABLE  5 

DEFINITIONS  OF  PARAMETERS  TRANSFERRED  TO  USTORE 


PARAMETER  DEFINITIONS 

I PLOT  THE  CURRENT  PLOT  NUMBER 

I RUN  THE  CURRENT  RUN  NUMBER 

IREC  THE  CURRENT  TAPE4  RECORD  NUMBER 

KREC  THE  DIMENSION  OF  THE  ARRAYS  TIME-'  AND  PI  THROUGH  P-NPS. 

KREC-  IS  THE  NUMBER  OF  RECORDS  (DATA  POINTS)  IN 
EACH  SOFE  RUN. 

NCRV  THE  NUMBER  OF  CURVES  TO  BE  DRAWN  FROM  DATA  STORED  IN 

THE  P  ARRAYS.  THE  ARRAYS  PI  THROUGH  P-NCRV  WILL  BE 
PLOTTED  AS  CURVES  ONE  THROUGH  NCRV.  0<NGRV<6  . 

NPS  THE  NUMBER  OF  P  ARRAYS  ALLOCATED  FOR  USER  PLOTTING 

AND  SCRATCH  STORAGE  OF  DATA.  NPS  MUST  BE  GREATER  THAN 
OR  EQUAL  TO  NCRV  AND  LESS  THAN  OR  EQUAL  TO  10.  NOTE 
THAT  ALTHOUGH  THE  ARRAYS  PI  THROUGH  P10  ARE  PASSED 
TO  THIS  ROUTINE  AS  FORMAL  PARAMETERS,  THERE  IS  NO 
STORAGE  ALLOCATED  TO  THE  ARRAYS  P-(NF'S+1)  THROUGH  P10. 

NS  THE  NUMBER  OF  ELEMENTS  IN  THE  TRUTH  STATE  VECTOR 

NF  THE  NUMBER  OF  ELEMENTS  IN  THE  FILTER  STATE  VECTOR, 

ALSO  THE  NUMBER  OF  ELEMENTS  IN  VF. 

XS  TRUTH  STATE  VECTOR 

XF  FILTER  STATE  VECTOR 

VF  STATE  VARIANCES  FROM  DIAGONAL  OF  THE  COVARIANCE  MATRIX  PF 

TIME  TIME  ARRAY,  DIMENSIONED  KREC.  ALL  VALUES  UP  TO  AND 

INCLUDING  TIME(IREC)  HAVE  BEEN  DEFINED  WHEN  USTORE 
IS  CALLED. 

P1-P10  LINEAR  ARRAYS  USED  TO  STORE  DATA  FOR  PLOTTING  AND  SCRATCH 
PURPOSES.  CURVES  WILL  BE  DRAWN  FROM  VALUES  IN  PI  THROUGH 
P-NCRV.  ARRAYS  PI  THROUGH  P-NPS  HAVE  BEEN  ALLOCATED  KREC 
WORDS  OF  STORAGE.  ARRAYS  P-(NPS+1>  THROUGH  P10  HAVE  NO 
STORAGE  ALLOCATED  TO  THEM  AND  THUS  SHOULD  NOT  BE  USED. 

THE  P  ARRAYS  FOR  WHICH  STORAGE  IS  ALLOCATED  ARE 
INITIALIZED  TO  ZERO  BY  SOFEPL. 
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TABLE  6 

DEFINITIONS  OF  PARAMETERS  TRANSFERRED  TO  USETPL 


PARAMETER  DEFINITIONS 

I PLOT  THE  CURRENT  PLOT  NUMBER 

RUNCNT  THE  NUMBER  OF  RUNS  PROCESSED  TO  CREATE  THIS  PLOT,  AS 

SPECIFIED  ON  THE  ’RUN  SELECTION-'  CARD.  RUNCNT  IS  USEFUL 
FOR  COMPUTING  AVERAGES  AND  STANDARD  DEVIATIONS. 

KREC  THE  DIMENSION  OF  THE  ARRAYS  'TIME-'  AND  PI  THROUGH  P-NPS. 

KREC  IS  THE  NUMBER  OF  RECORDS  (DATA  POINTS)  IN 
EACH  SOFE  RUN.  . 

NCRV  THE  NUMBER  OF  CURVES  TO  BE  DRAWN  FROM  DATA  STORED  IN 

THE  P  ARRAYS.  THE  ARRAYS  PI  THROUGH  P-NCRV  WILL  BE 
PLOTTED  AS  CURVES  ONE  THROUGH  NCRV.  0CNCRVC6  . 

NFS  THE  NUMBER  OF  P  ARRAYS  ALLOCATED  FOR  USER  PLOTTING 

AND  SCRATCH  STORAGE  OF  DATA.  NPS  MUST  BE  GREATER  THAN 
OR  EQUAL  TO  NCRV  AND  LESS  THAN  OR  EQUAL  TO  10.  NOTE 
THAT  ALTHOUGH  THE  ARRAYS  PI  THROUGH  PlO  ARE  PASSED 
TO  THIS  ROUTINE  AS  FORMAL  PARAMETERS,  THERE  IS  NO 
STORAGE  ALLOCATED  TO  THE  ARRAYS  P-(NPS+1)  THROUGH  PlO. 

TIME  TIME  ARRAY,  DIMENSIONED  KREC. 

PI -PlO  LINEAR  ARRAYS  USED  TO  STORE  DATA  FOR  PLOTTING  AND  SCRATCH 
PURPOSES.  CURVES  WILL  BE  DRAWN  FROM  VALUES  IN  Pi  THROUGH 
P-NCRV.  ARRAYS  PI  THROUGH  P-NPS  HAVE  BEEN  ALLOCATED  KREC 
WORDS  OF  STORAGE.  ARRAYS  P- ( NPS+ 1 )  THROUGH  PlO  HAVE  NO 
STORAGE  ALLOCATED  TO  THEM  AND  THUS  SHOULD  NOT  BE  USED. 

THE  P  ARRAYS  FOR  WHICH  STORAGE  IS  ALLOCATED  ARE 
INITIALIZED  TO  ZERO  BY  SOFEPL. 


be  discussed  briefly  in  Section  5  for  an  example  problem.  The 
Section  5  example  illustrates  the  essentials  of  acquiring  the 
appropriate  files,  operating  on  them  to  create  an  executable  load 
module,  and  then  running  the  program  to  produce  a  PLFILE.  The 
example  JCL  is  easily  modified  to  accomodate  a  user's  special 
needs,  such  as  substituting  magnetic  tape  for  disk  input  of 
TAPE4.  The  procedures  for  such  modifications  are  discussed  in 
the  appropriate  CDC  and/or  ASD  Computer  Center  manuals  (eg.  [4] 
and  [7])  and  are  not  covered  here  except  to  note  the  following. 
Communications  into  and  out  of  SOFEPL  are  accomplished  via  the 
data  files  listed  on  the  PROGRAM  card,  the  first  card  in  the 
SOFEPL  source  deck.  This  card  is  a  special  CDC  construction  that 
reads  as  follows: 


PROGRAM  SOFEPL  (TAPE5  =  64/80,  TAPE4,  OUTPUT,  TAPE6  =  OUTPUT, 
PLFILE  =  0) 

TAPES  is  the  plot  specification  file,  TAPE4  the  SOFE/SOFEPL 
interface  file,  and  PLFILE  the  device  independent  plot  file. 
OUTPUT  and  TAPE6  are  included  on  the  PROGRAM  card  to  accomodate 
printed  output  requirements. 


SECTION  5 
EXAMPLE  RUN 


This  section  demonstrates  SOFEPL's  capabilities  with  an 
example.  Seven  plots  are  created  for  the  SOFE  satellite  orbit 
example  discussed  in  [1] .  Understanding  the  example  presented 
here  does  not  require  a  detailed  knowledge  of  the  SOFE  example, 
but  it  will  be  helpful  to  recall  these  facts.  The  orbit  problem 
is  expressed  in  polar  coordinates  with  truth  and  filter  state 
vectors  of  dimension  four.  The  state  vector  quantities  are,  in 
order,  range,  range  rate,  angle,  and  angle  rate.  These  facts  are 
reflected  in  plot  titles  and  in  SOFEPL  control  parameter  choices. 
SOFEPL's  printed  output  and  all  seven  plots  appear  in  Appendix  A. 

5 . 1  Example  Input 

Figure  9  shows  the  JCL  commands,  the  Update  directives, 
and  the  FORTRAN  code  required  to  run  SOFEPL  for  the  example 
problem.  The  first  line  of  JCL  in  Figure  9  is  the  job  card  that 
adheres  to  a  format  established  in  [4] .  The  next  three  lines  are 
comments  that  are  echoed  in  the  dayfile.  In  the  fifth  line,  the 
SOFEPL  source  code  is  ATTACHed  (brought  to  the  job  from  disk 
storage  as  a  local  file  named  OLDPL) .  The  SOFEPL  source  code  (on 
the  file  OLDPL)  is  maintained  in  a  compressed  format  using  CDC ' s 
Update  utility  [5] .  The  sixth  command  is  an  Update  control  state¬ 
ment  that  expands  OLDPL  according  to  the  Update  directives 
placed  after  the  *EOR  in  Figure  9.  This  expansion  creates  a 
file  named  COMPILE  that  is  amenable  to  FORTRAN  compilation.  The 
COMPILE  file  so  created  contains  all  the  subroutines  of  SOFEPL, 
with  USTORE  and  USETPL  modified  by  the  FORTRAN  code  following  the 
two  * INSERT. . .  directives.  These  modifications  are  needed  for 
the  desired  type-18  plot. 
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SHM.  T35.  1060.  CM1 15000.  problem  no.  .  n4,me»  phone,  etc. 
COMMENT. 

COMMENT.  *♦  SAVED  IN  SOFEPLJCL, ID«SHM. CY-999, SN*AFAL  ** 
COMMENT. 

ATTACH. OLDPL. SOFEPL.  ID=SHM,  MR«1 ,  CY*999,  SN*AFAL. 

UPDATE , F , C=COMP I LE . 

FTN,L=0. I*COMPILE. 

RETURN , OLDPL . COMP I LE . 

ATTACH. DISSPLA, DISSPLA, I D*L I BR ARY »  SN**ASD »  CY=999 . 

L I BRARY , D I SSPL A . 

ATTACH. TAPE4, SOFEORB I TPLOTT APE , ID=SHM. CY*1 , SN*AFAL. 
ATTACH, DATA, SOFEPLDAT A, ID«SHM, MR=1 , CY*999, SN»AFAL. 
REQUEST , PLF I LE , *PF . 

LOO < DATA) 

CATALOG , PLF I LE , PLF I LE , I D«SHM . 

♦EGR  <END  OF  RECORD  MARK) 

♦ID  EXAMPLE 
♦INSERT  UST0RE.62 
C 

C  CONVERT  RANGE  AND  ANGLE  COORDINATES 
C  TO  CARTESIAN  COORDINATES 

C 

Pi (IREC)  =  PKIREC)  +  XS<1)  *  C0S<XS<3>> 

P2 <  IREC )  =  P2(  IREC)  ♦  XS(1)  *  SIN(XSO)) 

P3( IREC)  ■  P3CIREC)  *  XF(1)  *  C0S(XF(3>) 

P4 ( IREC >  -  P4 ( IREC )  ♦  XF(1)  *  SIN(XF(3)) 

♦INSERT  USETPL.47 
C 

C  COMPUTE  AVERAGE  AT  EACH  POINT  OVER  ALL  RUNS 
C 

DO  10  I  *  l.KREC 

PI  ( I )  =  PHD  /  RUNCNT 
P2< I )  =  P2< I )  /  RUNCNT 
P3( I )  =  P3< I >  /  RUNCNT 
P4 ( I )  ■  P4(I)  /  RUNCNT 
10  CONTINUE 


\ 


JCL 


UPDATE 
Directives 
and  Code 
for  Plot 
Type  18 


Figure  9.  JCL,  UPDATE  Directives,  and  Fortran 
Code  for  the  Example  Run. 
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In  line  seven  of  the  JCL  in  Figure  9,  the  FORTRAN  compiler 
[61  produces  an  object  file  named  LGO  from  the  source  code  in 
COMPILE.  The  L=0  setting  suppresses  listing  of  the  code.  In  the 
eighth  JCL  line,  OLDPL  and  COMPILE  are  RETURNed  in  order  to  free 
disk  space  and  reduce  system  overhead.  In  lines  nine  and  ten, 
the  library  of  DISSPLA  routines  is  ATTACHed  and  formally  declared 
a  LIBRARY.  In  lines  11  and  12,  TAPE4  and  the  plot  specification 
file  (called  DATA)  are  ATTACHed  from  disk  locations  where  they 
have  previously  been  CATALOGed.  Line  13  requests  space  on  disk 
for  the  'permanent'  storage  of  PLFILE.  Line  14,  LGO ( DATA ) , 
causes  the  computer  to  load  and  execute  SOFEPL  using  the  files 
DATA  and  TAPE4  for  input  and  creating  PLFILE  as  output.  Line  15 
CATALOGS  PLFILE  for  subsequent  input  to  a  DISSPLA  postprocessor. 

It  is  beyond  the  scope  of  this  document  to  explain  further 
the  JCL  and  Update  commands  used  in  Figure  9.  Users  needing  more 
information  are  referred  to  the  references  cited  above  and  to 
Reference  7  which  covers  the  CDC  operating  system. 

As  noted  above,  the  SOFEPL  plot  specification  file  is  the 
local  file  named  DATA.  Figure  10  shows  this  file,  consisting  of 
a  title  and  seven  plot  sets.  Blank  lines  are  inserted  between 
plot  sets  to  enhance  readability.  Such  lines  are  not  required, 
but  they  do  no  harm  when  used  as  shown.  Discussion  of  the  plots 
themselves  is  in  Subsection  5.3. 

5 . 2  Example  Printed  Output 

Appendix  A  contains  11  pages  of  printed  output  from  the 
example  run.  The  first  two  pages  are  system  responses  to  the 
control  statement  UPDATE,  F,  C=COMFILE.  These  two  pages  show  the 
code  inserted  in  USTORE  and  USETPL,  and  summarize  the  deck  struc¬ 
ture  of  the  OLDPL  file.  The  third  page  is  the  introduction  page 
that  marks  the  beginning  of  SOFEPL  execution.  This  page  gives 
the  title,  the  run  date  and  time,  and  information  about  the  SOFE 
run  taken  from  the  TAPE4  header. 

Pages  four  through  six  of  the  printed  output,  starting  with 
the  header  SOFEPL  CARD  DATA  INPUT  AND  STORAGE  ALLOCATION,  echo 
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SOFEPL  EXAMPLE  BASED  ON  SATELLITE  ORBIT  PROBLEM 

3,2.2/ 

1-L 

TRUTH  STATE  2  AND  FILTER  STATE  2  -  PLOT  TYPE  3* 
TIME  (TIME  UNITS)* 

RANGE  RATE  (LENGTH  PER  UNIT  TIME)* 


6,2.2/ 

1-L 

RANGE  RATE  ERROR.  AUG.  AND  ACTUAL  STD.  DEV. ,  PLOT  TYPE  6* 
TIME  (TIME  UNITS)* 

RANGE  RATE  ERROR  (LENGTH  /  UNIT  TIME)* 

12,1,1,1/ 

1-10, 20-30 , 40-50 

STD.  DEV.  OF  ACTUAL  AND  PREDICTED  EST.  ERROR  -  TYPE  12  PLOT 
TIME  (TIME  UNITS)* 

RANGE  STANDARD  DEVIATION  (LENGTH)* 

14,4,4,4,  ,  ,  ,  ,  57 ■  <i/ 

1-L 

COMPOSITE  ERROR  PLOT,  SAMPLE  OF  PLOT  TYPE  14* 

TIME  (TIME  UNITS)* 

ANGLE  RATE  ERROR  (DEGREES  PER  UNIT  TIME)* 

16, ,,,1,1,1/ 

1-L 

RANGE  RESIDUAL,  AVG.  AND  PREDICTED  STD.  DEV.,  PLOT  TYPE  16* 
TIME  (TIME  UNITS) 

RESIDUAL  (LENGTH  UNITS) 

18. . 4.4/ 

41-50 

TRAJECTORY  IN  CARTESIAN  COORDINATES  -  PLOT  TYPE  18* 

TIME  (TIME  UNITS)* 

X  AND  Y  POSITION,  TRUE  AND  ESTIMATED  (LENGTH  UNITS)* 

3.2.2. .  . . . 10. .3. .0.0.25. »-. 25, +. 25/ 

1-L 

SAME  AS 'PLOT  1  BUT  WITH  AXIS  SCALING  AND  TIME  WINDOWING* 
TIME  (TIME  UNITS  *  10.)* 

RANGE  RATE  (LENGTH  PER  UNIT  TIME  *  3. )* 


Figure  10.  Plot  Specification  File  for  the  Example  Run. 
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the  input  data  for  each  plot  and  show  the  starting  location  and 
memory  required  to  produce  each  plot.  Error  messages  may  appear 
in  this  section.  S0F5PL  attempts  to  recover  gracefully  from 
errors,  usually  by  skipping  plots  in  which  input  errors  are 
detected . 

Pages  seven  through  ten,  starting  with  the  header  PLOTTING 
EXECUTION,  are  output  after  the  data  for  each  plot  is  generated. 
These  pages  contain  a  summary  of  important  data  for  each  plot, 
some  of  it  a  repeat  of  previously  echoed  input  data  placed  here 
for  conveninece,  and  some  of  it  new  data  available  only  when  all 
curves  of  a  plot  are  generated.  In  this  latter  category  are  the 
actual  number  of  points  in  each  curve,  min/max  values  on  each 
axis,  and  the  average  of  the  y-data  for  each  curve.  This  section 
of  output  also  contains  a  listing  of  the  t,  y  coordinates  for 
each  curve  if  IPNTSW=1  {see  plot  number  5) .  The  last  page  of 
printed  output  is  the  job  dayfile  in  which  the  JCL  commands  of 
Figure  9  appear  along  with  the  system's  responses. 


5 . 3  Example  Plots 

The  seven  plots  produced  by  the  example  run  are  shown  as 
the  last  seven  pages  of  Appendix  A.  The  first  five  plots 
demonstrate  a  representative  cross-section  of  available  SOFEPL 
plots,  while  the  last  two  focus  more  on  SOFEPL* s  special 
capabilities.  Between  the  instructions  and  explanations  in 
Sections  2,  4,  and  5,  the  linkage  from  Figure  10  through  the 
printed  output  to  each  plot  should  be  clear.  Note  the  following 
about  the  first  group  of  five  plots:  in  Plot  3,  a  selected  set 
of  32  of  the  available  50  runs  are  processed;  in  Plot  4,  angles 
are  converted  from  radians  to  degrees  using  57.3  for  SFY;  in  Plot 
5,  both  axis  lables  are  oft -center  because  terminating  $  signs 
were  not  used  on  input. 

Plot  6  illustrates  plot  type  18,  the  user-defined  case.  In 
this  example,  the  orbiting  satellite's  r-0  coordinates  are 
expressed  in  x-y  Cartesian  form.  To  accomplish  this,  the  code 
shown  in  Figure  9  is  inserted  in  USTORE  and  USETPL.  This  code 


forms  the  x-  and  y-positions,  both  estimated  and  true,  and  aver¬ 
ages  them  over  the  last  10  of  the  available  50  runs.  Note  that 
each  curve  in  Plot  6  has  a  different  line  pattern  to  help 
distinguish  it  from  its  neighbors. 

Plot  7,  the  last  plot,  demonstrates  the  use  of  axis  limits 
and  scale  factors.  In  this  plot  these  parameters  have  been  cho¬ 
sen  so  as  to  "zoom  in"  on  the  data.  Plot  7  should  be  compared 
with  Plot  1  to  see  the  zooming,  or  windowing,  effect.  Both  plots 
are  examples  of  plot  type  three  and  both  use  the  data  from  the 
same  runs.  The  scale  factors,  axis  limits,  and  titles  are  the 
only  differences.  Note  that  the  axis  limits  are  in  the  units  of 
the  data  after  the  scale  factor  was  applied. 

Different  line  patterns  are  used  in  plots  with  multiple 

curves  to  distinguish  which  curve  corresponds  to  what  quantity. 

In  plot  type  18  each  curve  has  a  unique  line  pattern.  In  plot 

types  1  through  17,  these  rules  apply.  When  only  one  curve  is 

drawn,  its  line  is  solid.  When  two  curves  are  drawn,  the  line 

pattern  is  solid  for  the  truth  quantity  and  broken  for  the  filter 

quantity.  When  three  curves  are  drawn,  two  of  them  are  always 

symmetric,  either  about  y-o  or  about  the  third  curve;  in  this 

case  the  symmetric  pair  is  broken  and  the  solitary  curve  is 

solid.  There  are  no  four-curve  plots.  In  the  single  five-curve 

plot,  three  line  patterns  are  used:  E.  .  is  solid;  E. .  +  SE..  are 

.  ■  i]  13  —  13 

broken  of  one  pattern;  and  +VVf .  are  broken  of  another  pattern. 
In  this  last  case  differences  in  symmetry  help  to  distinguish 
which  curves  are  which. 


\ 


APPENDIX  A 
EXAMPLE  RUN  OUTPUT 


o 

<L 


•r 

N 


N 

•N 


oo 

04 

UO 

I 

«* 


Ui 

H 

<1 

Q 

LL 

D 


»N  >0  *  JO  >0  N  00 


o-nNro^sr>-opv 


Ui  UJ  UI  UJ  UJ  UJ  UJ  UJ  UJ 
Q-OLii-Q_a,u_U.d.a. 

zee  eee  zee 

xxxxxxxxx 
UllilUl  UiUJUJ  UiUJUJ 


JJ  UJ  LU  LU  UJ  UJ  UJ  LU 
-J-l-l-J-l-l-J-J 

zzzclzzzz 
eez  eez  ee 
^^<x<x<x<i*x<£ 

xxxxxxxx 
UJ  UJ  UJ  UJ  UJ  UJ  LU  UJ 


•— 

/%  />  « 

CO 

<0 

*■>  a  « 

z 

z 

A  A  A  A 

00 

oo  o  oo 

D 

D 

OO  00  00  00 

Z 

iZ 

0 

UJ  U.  LL 

-J 

_J 

x  X  x  X 

-J 

_l 

•A 

z  0  z 

<1 

<1 

0  Z  0  Z 

0 

o 

*-•  o  *-« 

CO 

o  *-»  O  *-• 

LU 

UJ 

o 

CO  w  «0 

z 

(Z 

UJ 

~i 

h- 

LU 

UJ 

K 

Z 

<1 

* 

*  *  * 

> 

> 

<1 

*  *  *  * 

i- 

Z 

o 

a 

z 

<3 

i-* 

^  ^ 

A  A  A  M 

X 

a 

•H 

^  H  ^ 

H- 

Q 

Ui 

Z 

z 

«-  H  f-  J- 

Z  Hh 

H  K  iZ 

V  s„  V  V 

o 

iO 

0  LL  Ll 

►4 

z  z  z  z 

-  z  z 

z  z  o 

o 

X 

XXX 

o 

g  u  o  u 

2 

o  oo 

O  O  o 

•J 

LL 

zzzz 

i3 

CL  Z  Z 

z  z  o 

CO 

+ 

+  +  + 

3DD3 

Z 

Z>  3 

3D  vO 

+  +  4-  4* 

t- 

UJ  UJ 

I 

Z  Z  Z  Z 

«X 

X  £Z  a: 

£Z  iZ  UJ  UJ 

z 

-J  H- 

O 

a 

o 

-J  h- 

A  A  A  A 

IU 

s  s 

o 

O  w  u 

<r 

S  N  N  V 

<1  ^  v 

N  \  o  C 

a 

Z  Z 

UJ 

LU  LU  LU 

UJ 

-j 

UJ 

z  z 

UJ  UJ  UJ  0 

** 

<1  — 

a: 

z  z  a: 

o  ~  ~  ~  ~ 

o 

Q 

ZZZZ 

Q 

M  M  H 

K 

LU  M  H  M  M 

Z 

K  UJ  «  i-i 

»t  M  Q 

a  Z 

<1 

z  w  ~  ^  ~ 

h* 

<1  QC  — 

—  ^  Q  tZ 

w  w  w  ^ 

4  2 

<N  0  «}■ 

¥  •*  <N  00  * 

Z 

iC  **  04 

oo  ^  z  5 

-•  in  oo  «r 

a  o 

a. 

a.  z  a. 

UJ 

*  a.  z  a.  a. 

o 

uj  •  a.  a. 

a.  a.  3  o 

ll  rt 

CM 

o 

n 

a 

^  UJ 

u 

o  ^ 

UJ  o 

>0 

UJ 

a 

h  a  «  <r 

<i 

«  «  8  8  D 

<x  us 

u  n  D  Ui 

II  II  II  U 

• 

2  S 

• 

iZ 

H  Z 

V 

x  n 

Z  3  2 

UJ  LU 

Z  <c 

UJ 

^  ^  A  A  M 

UJ  — 

~  ~  z  <c 

A  A  A  A 

J  IT 

<r  ~ 

o 

<J  O  O  0. 

> 

(0 

>  M  «  —t 

a.  o 

£  CO 

UJ 

LU  UJ  UJ  >- 

«x 

A  - - -  Z 

z 

w  Z  ITu, 

E  K 

UJ 

L L 

Z  Z  Z  UJ 

o  c*  oo  <r  o 

5 

O  <N 

oo  o  uj 

ZZZZ 

<x  0 

H  H 

•H  **  M  (J) 

UJ 

-•  o.  a.  a.  a.  o 

UJ  —  Cl  a. 

0-  Q_  u  k. 

x  D 

Z  Z 

r  -  -  3 

H 

h* 

V- 

QC  QC 

V  w  W  L^ 

UJ 

LU  <X 

«-4 

<N«^ 

D 

o 

<i 

3  O 

UJ  <x 

-(N  CO  ♦ 

h 

>  O 

a. 

a.  a.  a.  t- 

Of 

Q 

u 

a.  a 

>  U 

ZZZZ 

-J 

H-  CC 

z 

Z 

£ 

E 

z 

a. 

Z  Ul 

o  o 

3 

a. 

O 

o  o 

M 

UJ  0 

CJ  J— 

0 

o 

M 

(j 

U  K 

W 

Q  Z 

z 

a 

o 

*X  P“» 

*■4 

O 

Q 

o 

a 

*  *  o  o  o 

u 

* 

c 

c 

c 

** 

£ 

DUO 

—  u  o  u 

<J 

-j 

UJ 

J  J  J  J  J  J 

-J  -1  -J  UJ  UJ  UJ  ui  UJ  UJ  UJ  UJ 

a 

<x 

s  * 
s  * 

* 

* 

U0.IUQ. 

Kh(-KhK 

Sggggggg 

-1 

'v  * 

♦ 

UJ  UJ  LU  m  14!  UJ 

UIUJUJ  t-HKH-Kk-KK 

•>.  * 

* 

0  </5  Cfi  (0  CO 

CO  0  U)  0  0  CO 

3 

^  * 

* 

333335 

333  333 

3  3  5  3  3 

42 


*#***•*#**«#**#**•«*«*-»•*  SOFEFL  CAKD  DATA  INPUT  AND  STORAGE  ALLOCATION  «*««««««***#**«»*«*««««*#««*#«*#««««* 


PLOT  4  INPUT  INFORMATION: 


*************************** **********  PLOTT INU  EXECUTION  *************************************************** 


I 


UJ 

CL 

> 


O 

-1 

CL 


3> 

UJ 

Q 


Q 

H 


Lli 

cc 

LL 


z 

<r 


o 

> 

«i 


<L 

3 

a 


-S3 

a.  q. 

ZU 


u) 

UJ 

CC 


UJ 

3 

z 

sn <t 
ic 


<n 


«n  i 


vhL 

u>  oc  a: 

uj  o  q 

(C  U>  *f> 


%'i 

3U1 
-J 
UJ  w 

z 


vOO 

o>> 

>0  CO  >  «  « 

-  «  ♦  I 


OOt-C 
0  0-3 
O  O  Q 


o 

in 


^  Z 

**  Ui 
H  iX 


3S 

U.  K 


_  38 

ui  «  _i  h-  a 

<^8  5 S 

y  ui  «  j  3  r 


u 

> 


S  5  -*  »N  o) 


UJ  UJ  UJ  UJ 

H  »-  (-*>>> 

33  35555 

a.  a.  o.  z  u  u  u 


(/)  U  >A  <1  Z  -» 

w  «  «j  **  z 
x  x  iX  3 
<t  3  <1  0)  CL  (K 

I  H  |  M 

SSSSStf 

£>?;ss 


00 

04 

a> 

CN 

O 

>0 

«T 

-0 

04 

* 

«r 

O' 

fv 

O 

04 

r- 

O 

'O 

04 

<Y> 

0) 

00 

' 0 

*■ 

to 

04 

■4) 

>0 

00 

N 

U3 

04 

O 

<*k 

0 

0 

*■ 

C4 

04 

04 

M 

04 

04 

M 

04 

04 

Oil  III  I 


I  I  I 


or'-<<NC£>ovoa>-*.o»* 
N4N4<tO)  ■<  »  »  s  O 
NNO-4«Ni»niM(B 

>M)»iN«00000 

^♦iNNCMN  NNNNN 


% 

u 


<N 

O 

I 

UJ 

CO  00  CM  CM 
O  O  s)  •O  S  'O  -0 
OOMNC9  'O-O 
O  G*0*0fh't'+ 
o  O  N  N  ^  ^ 

in  •  «*  *  a>  C4  «n 

•  in . 

-  I  C4  I 


04  04  *  CM  04 

ooooo 

I  I  I  I  I 
UJ  UJ  UJ  UJ  UJ 
04  i>  •*  ^ 
— <  GO  — •  -0  CN  O* 
04  m  (N  -•  ♦ 
UJ  *  *  C4  *  © 
>  ^  10  O'  o 

a:  ^  O'  «r  in  o 
3  ■  •  ■  •  • 
O  0)  ^  i>  (0  04 


C4  *N  CO 

ooo 

I  i  I 
UJ  UI  UJ 
•OMN 
«> 

«r  oo  ■** 

rN  c«)  — 


CO  >0  ( 


04  CM 
O  O 
I  I 
UJ  UJ 
O'  * 
00  CO 
*  00 
o  00 
00  o 
» <r 


«o  co  in  *4)  'O 


UJ  UJ  UJ 

r  2  r 


UJ  UJ  UJ 


^S^^I>>ooo 

4  —  —  >  >  > 

WhK>>4 44 

a. 

slllhv* 

£s2ridud 


o  o  * 

UJ  o  o  < 

£OQ< 
wOOl 

h  in  • 

•  ^  W  04  04  Wflf  *  tn 


>000000 

>000000 

JS  8SS  gi 


2-hNk)»|(14n00».O 


2 


50 


I 


PLOT  NUMBER  »  6 

PLOT  TITLE  s  TRAJECTORY  IN  CARTESIAN  COORDINATES  -  PLOT  TYPE  IS* 


TRUTH  STATE  2  AND  FILTER  STATE  2  -  PLOT  TYPE  3 


SO'O 


(TIME  UNITS) 


nut  J7  flue,  laai  jab-9ranr  ,  wnrs/nao  biwKfl  vw  r.j 


0 


2.0  2.5  3.0  3.5  4.0  4.5 

TIME  (TIME  UNITS) 


«3A  WTdSSlO 


56 


(TIME  UNITS) 


RANGE  RESIDUAL,  AVG.  AND  PREDICTED  STD.  DE 


I  ME 


TRAJECTORY  IN  CARTESIAN  COORDINATES  -  PLOT  TYPE  18 


(TIME  UNITS) 


APPENDIX  B 

DISSPLA  POSTPROCESSORS 


The  DISSPLA  plotting  routines  referenced  within  SOFEPL  out¬ 
put  a  stream  of  compressed,  device  independent  plot  commands  on 
the  file  PLFILE.  It  is  the  job  of  the  DISSPLA  postprocessors  to 
transform  these  generalized  plot  commands  into  plot  instructions 
for  a  particular  plotting  device.  Thus  the  plot  data  from  one 
SOFEPL  run  stored  on  PLFILE  can  be  plotted  on  many  different 
plotting  devices  simply  by  running  the  appropriate  DISSPLA  post¬ 
processor.  Two  widely  used  postprocessors,  ONLINE  for  an  on¬ 
line  Calcomp  plotter  and  TEK4014  for  an  interactive  Tektronix 
4014  terminal,  are  reviewed  here. 

The  ONLINE  postprocessor  is  used  to  translate  the  DISSPLA 
plot  commands  into  commands  useable  by  an  on-line  Calcomp 
plotter.  This  postprocessor  can  be  executed  from  any  terminal  or 
batch  job.  ONLINE  uses  the  local  file  PLFILE  as  input  and 
creates  an  output  file  called  PLOT  containing  the  Calcomp  plot 
commands.  This  file  is  then  routed  to  the  on-line  plot  device. 
The  following  commands  can  be  used  to  execute  ONLINE. 

ATTACH,  DISSPLA,  ID=LIBRARY,  SN=ASD. 

LIBRARY,  DISSPLA. 

ATTACH,  PLFILE,  ID=SHM. 

ONLINE. 

ROUTE,  PLOT,  DC=PU ,  FID=SHM,  TID=AG . 

In  this  use  the  plot  file  produced  by  the  ONLINE  postprocessor 
was  routed  to  terminal  AG  to  be  plotted  on  the  AG  (Building  22) 
online  Calcomp  drum  plotter.  Upon  typing  the  ONLINE  command,  the 
user  will  be  prompted  to  enter  DISSPLA  postprocessor  directives. 
These  directives  and  their  effects  are  explained  later  in  this 
section. 

The  TEK4014  postprocessor  is  used  to  draw  the  plots 
directly  on  a  Tektronix  4014  Computer  Display  Terminal.  TEK4014 
is  an  interactive  processor  executed  while  the  user  is  logged  in 


on  the  Tektronix  terminal.  The  commands  used  to  execute  TEK4014 


are 


ATTACH,  DISSPLA,  ID=LIBRARY,  SN=ASD. 

LIBRARY,  DISSPLA. 

ATTACH, PLFILE,  ID=SHM. 

TEK4014 . 

The  postprocessor  then  prompts  the  user  for  the  needed  direc¬ 
tives.  After  the  directives  have  been  entered  the  postprocessor 
draws  the  plots  on  the  screen  one  at  a  time,  pausing  after  each 
plot  for  inspection  or  hard  copying  by  the  user. 

Plot  directives  instruct  the  DISSPLA  postprocessor  to 
selectively  plot,  scale,  window,  and  orient  the  ensuing  plots  on 
the  actual  plotting  surface.  The  plot  directives  are  the  same 
for  all  the  different  DISSPLA  postprocessors.  The  two  most  use¬ 
ful  directives  are  the  DRAW  and  the  MODIFY  commands.  The  DRAW 
command  specifies  which  plots  are  to  be  drawn  and  has  the  form 

DRAW=list 

where  "list"  is  the  list  of  plots  or  series  of  plots  to  be  drawn. 
For  example: 


DR.\W=1 ,3,5-7 , 10-13 , 16* 

would  araw  plots  1,3,5,6,7,10,11,12,13,16.  The  MODIFY  command 
controls  changes  to  the  plots  to  be  drawn  and  must  follow  the 
DRAW  command.  The  MODIFY  command  has  the  form 

MODIFY=list(  sped* spec 2*  . .  .  *specn) 

where  "list"  is  the  list  of  plots  to  be  modified  and  "spec"  is  a 
command  defining  a  change  to  be  made  to  the  plots.  Several 
specifications  (spec)  are  available  to  change  the  size,  position 
or  portion  of  the  plot  to  be  drawn.  Two  such  specifications  are 


SCALE  and  CORNER.  A  specification  of  SCALE=.75  acts  to  reduce 
the  plot  to  75  percent  of  its  original  size.  The  specification 
CORNER=1.5,1.5  would  move  the  lower  left  corner  up  and  to  the 
right  1.5  inches.  All  postprocessors  continue  to  read  directives 
until  a  $,  signifying  end  of  directives,  is  encountered.  A  set  of 
directives  to  draw  the  first  seven  plots  on  PLFILE,  reduce  their 
size,  and  recenter  them  is 

DRAW=7* 

MOD IF Y=  7 ( SC ALE= .75)* 

M0DIFY=7 (C0RNER=1 .5,1.5)  $ 

Note  that  the  "*"  is  necessary  to  separate  the  individual 
commands. 

A  detailed  discussion  of  all  DISSPLA  postprocessor  commands 
is  found  in  Reference  9. 
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APPENDIX  C 

MERGING  INTERFACE  FILES 


It  is  occasionally  required  that  two  or  more  'similar' 
interface  files  (TAPE4s)  be  merged  to  form  a  single  interface 
file  so  that  SOFEPL  plots  can  be  generated  from  a  composite  data 
set.  Two  'similar'  interface  files  would  exist,  for  example,  if 
two  SOFE  studies  were  made  of  the  same  problem,  one  study  of 
N i  runs  from  ISEED^ ,  and  a  second  study  of  N2  runs  from  ISEED2 • 
Merging  all  runs  into  one  master  file  of  N^  +  N2  runs  produces  a 
larger  ensemble  from  which  to  draw  conclusions. 

In  a  second  case,  'similar*  interface  files  are  also 
generated  when  one  SOFE  study  is  the  time-continuation  of  a  pre¬ 
vious  study.  In  this  case,  each  study  contains  the  same  number 
of  runs  but  the  later  study  takes  up  in  time  where  the  first 
ends.  SOFE  provides  this  capability  through  its  ICONT  parame,  r. 
Again,  benefits  could  derive  from  merging  the  TAPE4  outputs  of 
these  two  studies  into  one  file  of  longer  runs. 

This  appendix  discusses  two  programs  for  accomplishing  file 
merging.  The  first  is  named  AFFIXR  and  deals  with  merging  two 
similar  files  containing  runs  over  the  same  time  interval.  The 
second  is  named  AFFIXT  and  deals  with  merging  two  similar  files 
of  equal  numbers  of  runs  where  the  second  file  is  the  time- 
continuation  of  the  first. 

Breaking  a  large  study  into  two  or  more  substudies  may  toe 
advisable  for  any  of  these  reasons:  computer  turn-around  times 
for  big  runs  are  prohibitively  long;  a  short  preliminary  evalu¬ 
ation  of  performance  is  needed  before  a  lot  of  computer  time  is 
spent  and  possibly  wasted;  the  analyst  is  unsure  of  the  number  of 
runs  (or  run  lengths)  required  to  reach  statistical  significance 
so  he  would  prefer  to  monitor  convergence  rates  as  new  runs  are 
added  gradually. 


■■  sc. 
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AFFIXR 


AFFIXR  (affix  runs)  merges  two  interface  files  that  meet 
the  following  'similarity'  criteria: 

•  NS,  NF,  M  same  in  both  files 

•  TO,  TMEASO,  TF  same  in  both  file.* 

•  DTMEAS ,  DTCCPL  same  in  both  files 

The  IPASS  value  (number  of  runs)  of  each  file  may  be  different. 

AFFIXR  is  designed  to  execute  from  an  interactive  terminal 
connected  to  the  ASD  CDC  computer  system.  The  formal  names  of 
AFFIXR 's  input  files  are  TAPE1  and  TAPE 2 ,  and  that  of  its  output 
file  is  TAPE3 .  If  we  assume  that  the  file  AFFIXR  is  a  loadable 
object  module,  the  program  can  be  run  by  entering 

AFFIXR,  TAPE1 ,  TAPE2 . 

AFFIXR  checks  both  input  TAPES  for  compatibility,  i.e.,  checks  to 
see  if  they  satisfy  the  similarity  criteria  listed  above.  If 
they  are  not  compatible,  AFFIXR  prints  diagnostic  messages  as  to 
why  they  cannot  be  merged.  If  they  are  compatible,  the  program 
creates  a  TAPE3  output.  TAPE3  is  a  valid  data  file  for  input  to 
SOFEPL. 

AFFIXR  performs  its  merge  operations  in  the  following 
order.  First,  TAPEs  1  and  2  are  rewound  and  their  header 
records  are  read  and  compared.  If  the  TAPES  are  compatible, 

AFFIXR  writes  a  new  header  to  TAPE 3  in  which  IPASS3  =  IPASS^  + 

IPASS2.  This  header  also  contains  the  data  listed  in  the  simi¬ 
larity  criteria  plus  the  run  date,  time  and  title  from  TAPE1 . 

Second,  all  of  the  runs  from  TAPE1  are  copied  to  TAPE3 .  Third, 
all  the  runs  from  TAPE 2  are  copied  to  TAPE3. 

Repeated  executions  of  AFFIXR  can  be  used  to  merge  any 
number  of  files.  Suppose  a  user  wants  to  merge  files  A,  B,  and 
C.  He  can  do  this,  with  the  final  output  on  E,  using  execution¬ 
time  replacement  of  file  names  as  follows: 

1 
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AFFIXR,  A,  B,  D 
AFFIXR,  D,  C,  E 


AFFIXT 

AFFIXT  (affix  time)  merges  two  interface  files  that  meet 
the  following  similarity  criteria: 

•  NS,  NF,  M  same  in  both  files 

o  TFX  <  T02 

•  IPASS  same  in  both  files 

The  second  criterion  states  that  the  two  studies  must  not  overlap 
in  time;  although  studies  may  be  separated  by  a  time  gap,  the 
usual  case  will  be  that  TF^  =  TO2. 

A  feature  of  AFFIXT  is  that  the  two  input  files  (TAPE1  and 
TAPE2  again)  do  not  have  to  match  in  DTMEAS,  DTCCPL  or  TMEASO. 
This  feature  is  useful  if  one  wants  to  merge  mission  segments  of 
different  types  to  obtain  a  composite  mission.  AFFIXT  is  exe¬ 
cuted  in  the  same  manner  as  AFFIXR  with  TAPE 3  again  being  the 
output  file. 

AFFIXR  and  AFFIXT  source  codes  can  be  attached  from  the  CDC 
disk  as  follows: 

ATTACH, AFFIXR, ID=SHM,SN=AFAL,CY=999 

ATTACH , AFFIXT , I D=SHM , SN=AFAL , CY=9  99 
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